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We investigate the quasiparticle spectrum near surfaces in a two-dimensional system with d-density- 
wave order within a mean-field theory. For Fermi surfaces with perfect nesting for the ordering wave 
vector (tt, n) of the d-density-wave, a zero energy bound state occurs at 110 surfaces, in close analogy 
with the known effect in d-wave superconducting states or graphite. When the shape of the Fermi 
surface is changed by doping, the bound state energy moves away from the Fermi level. Furthermore, 
away from half-filling we find inhomogeneous phases with domain walls of the d-density-wave order 
parameter. The domain walls also support low energy bound states. These phenomena might 
provide an experimental test for hidden d-density-wave order in the high-T c cuprates. 
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I. INTRODUCTION 

In the search for an explanation of the so-called 
pseudogap regime of the underdoped high-T c cuprates, 
Chakravarty, Laughlin, Morr, and Nayaktl recently pro- 
posed a novel type of order, known as d-density-wave 
order. This state was_discussed earlier by Schulza and in 
more detail by Nayaktl It breaks translational and time- 
reversal symmetry and is related to flux phases known 
from t-J-axid Hubbard models at infinite or large onsite 
repulsioELrtl Its order parameter, written with the elec- 
trons' creation and annihilation operators and ordering 
wavevector Q = (tt, tt), 



$ddw = 



k.s 



/(*) 



<i 



C k+Q,sl 



(1) 



is purely imaginary due to the sign change of the d x i_ y i 
wave form factor 



f(k) ~ cos k x — cos k y 



(2) 



under k — » k + Q. In real space it gives rise to circulat- 
ing in-plane currents around elementary plaquettes. The 
orientation of the current flow, i.e. the sign of the flux 
penetrating the plaquettes, is alternating from plaquette 
to plaquette. 

Since the staggered flux order parameter is of Ising 
type the d-density-wave can order in two spatial dimen- 
sions at nonzero temperatures. In one spatial dimension 
it can still lead to long range order in the ground state, 
as recently shown for the two-leg ladder by Fjaerestad 
and MarstonQ. 

Apart from the implications of the proposed type of 
order on observable bulk quantities Uke .-specific heat, 
spin susceptibilities and conductivitiesErQu, the d-wave 
structure of the order parameter opens another field of 
phenomena. These are surface effects, which have al- 
ready played a significant role in the experimental at- 
tempts to identify the superconducting pairing symme- 
tries in the high-T c cuprates and other unconventional 



superconductorsoEl Explicitly we mention the possi- 
bility of Andreev bound -stales at surfaces of uncon- 
ventional superconductorsErO with suitable orientation 
(and quality in experimental studies) at low or zero ener- 
gies. The theoretically appealing feature of these surface 
bound states is that they are a consequence of the sym- 
metry of the order parameter and depend less on the 
details of the (most of the time drastically simplified) 
models used to calculate them. Since the d-density-wave 
order described above also occurs with a non-trivial sym- 
metry, it appears to be interesting to search for surface 
effects in a d-density wave ordered state. Related to lat- 
ter phenomena are of course local deformations of the 
bulk state at impurities. In the context with d-density- 
wavc states, these were recently studied theoretically by 
Zhu et al.Ej, WangEj, and Mono, which all found low 
energy states in the vicinity of non-magnetic impurities. 

In this paper we analyze the possibility of surface 
bound states for electrons on a two-dimensional square 
lattice with nearest neighbor hopping in the presence of 
d-density-wave order. After a simplified analytical treat- 
ment showing the existence of surface bound states at 
[110] surfaces, we solve the mean-field equations for the d- 
density-wave amplitude self-consistently on a square lat- 
tice strip with finite width. The d-density-wave order is 
generated by a repulsive nearest neighbor density-density 
interaction. We find low energy surface states at the 
[110] surface of the lattice and study their dependence 
on the band filling. Further we analyze the d-density- 
wave state away from perfect (jr, 7r)-nesting at half fill- 
ing and find indications for domain walls of the order 
parameter forming between locally half-filled d-density- 
wave domains. These domain walls themselves can again 
host bound states with sub-gap energies. 

Of course a realistic description of the high-T c cuprates 
may require a more sophisticated theory, so that our 
simplified mean-field treatment might neglect a substan- 
tial part of the physics. Nevertheless we emphasize that 
the occurrence of surface bound states relies mainly on 
the symmetry of the order parameter. Thus our find- 
ings should be a robust feature remaining true in more 



detailed studies. Moreover the theoretical analysis of 
a large number of experiments in the framework of d- 
density-wave order by Chakravarty et al£l was performed 
on a comparable level of approximations. 



the fast oscillating and y-dependent parts from the wave 
functions by writing 



*A(B)(?) = fA(B), K (x')e i ™'e<«' 



(5) 



II. ANALYTICAL TREATMENT 

In order to understand the main features at the surface, 
let us consider electrons on a two-dimensional square lat- 
tice with hopping between nearest neighbors at half fill- 
ing. Then the free Fermi surface is the perfectly nested 
square. We assume that the system has d-density-wave 
order, giving rise to a term in the Hamiltonian, 

-ffddw = ^ iW (cos k x - cos k y ) 



C k c k+Q 



C- -Ct; 
k+Q k 



(3) 

after mean-field decoupling, with real Wq and Q = (-7T, it) 
(in units of the inverse lattice constant a~ l ) . In real space 
this term also couples nearest neighbor sites, analogous 
to the ordinary hopping term, so that the total Hamil- 
tonian has bipartite structure. In the following we will 
use this property and rewrite the Hamiltonian by dou- 
bling the unit cell and separating the wave function into 
components belonging to the A- and i?-sublattice, ^>A(r) 
and 4 , s(r), 







^k 



l 7fc 



* B (f) 



= E 



* A (r) 



(4) 



where e^ — 



Next we linearize the kinetic energy e^ around the line 
k' x = 7T and introduce p' x = k' x — it. Then, for fA(B).k' 
the Hamiltonian has the form 



H 
et 



where 



fA, K (x') 

!b, k {x') 



VF,x'(k'y)p' x 



= E 



fA,k> y (x>) 
fB,k>M) 



(6) 
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Here we replace p' x = —id x i. For the [110]-surface the 
outer most sites belong to one sublattice, say A. Then 
we assume that the wave function takes nonzero values 
only on this sublattice and vanishes on the other. This 
behavior of the bound state wave function is confirmed 
by the self-consistent solution described below. Thus, we 
impose the condition fs,k' (#') = and E — and ana- 
lyze the conditions to obtaip_.a zero energy bound state 
from the Hamiltonian abovecd The resulting equation is 



id x >fA,k>(x') = 2i 



VF,x 



sm 



fA,y(x') (7) 



-2t(cosk x + cos k y ) and 7^ = Wo(cosk x — 

cos k y ) and k = — zV. Obviously this formulation leads to 
the proper excitation spectrum in the uniform case which 
has two Fermi points on the boundary of the reduced 
Brillouin zone at k = K = (§, §) and K' = (§,-§). 
Close to these points the spectrum has the approximate 
form of a Weyl spectrum for a chiral particle with oppo- 
site chirality at K and K'. In this respect the spectrum 
of the d-density-wave state has |Obvious similarity with 
the spectrum of a graphite sheetEil. 

Now we consider the problem of a surface with nor- 
mal vector parallel to the [110]-direction which defines 
the coordinate x' — (x + y)J\/2 (x! > and the sur- 
face is located at x' = 0). Here we will not attempt 
a self-consistent solution of the mean-field problem, but 
keep the d-density-wave order parameter Wq constant for 
x 1 > 0. This is sufficient to show that this surface pos- _. 

sess surface bound states at zero energy. The transla- 
tion symmetry along the direction parallel to the surface 
(y' = (x — y)/\/2) is maintained so that the momentum 
k' y = (k x — ky)J\f2 is conserved. Let us consider elec- 
tronic states in the vicinity of the non-interacting Fermi 
line k' x — (k x + k y )/y/2 — -k (in units of inverse of the 
lattice constant a* = \[2a of the doubled unit cell) at 
the boundary of the reduced Brillouin zone. We extract 



which has an exponentially decaying bound state solution 
!a,w {x') = Ce~ KX . The existence condition is k' y > 
with k = 2W sm(k y /2) jv F ^ x and W > 0. As k' y ap- 
proaches the bound state extends deeper and deeper 
into the bulk. Zero energy bound states only exist for 
one sign of the momentum parallel to the surface, an in- 
dication of broken time reversal symmetry of the bulk 
phase. Furthermore, the bound state has non-bonding 
character as it occupies only one sublattice. Applying 
the time reversal operation or interchanging A- and B- 
sublattices reverses k' . A related situation occurs at so- 
called zigzag edges of graphite where the bound state 
spectrum has similar p r o p erties, but does not violate 
time reversal symmetrycSEl 
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FIG. 1. The [110] surface of the two-dimensional square 
lattice. Left: real space picture. The circles denote the 
lattice sites of A- and B-sublattice, and the differently 
shaded bonds between them indicate the sign structure of 
the d x i_ y 2 -density- wave amplitude 

fc-space view. 
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We may also include a small shift of the chemical po- 
tential /i (<C t). Within our approximative scheme this 
merely leads to a shift of the bound state energy E to 
— fi. A similar shift of-the spectrum was also obtained by 
Zhu et al.ta and Wangle for the impurity case. Below, we 
will find the same effect in the self-consistent treatment. 

Concluding this analysis we observe that the [110]- 
surface of a d-density-wave ordered state on a two- 
dimensional square lattice should support low-lying sur- 
face states. We note that at [100] surfaces no bound state 
can form according to a Bogoliubov-deGennes analysis. 
In the next section we confirm our qualitative results for 
the [110] surface with a self-consistent mean- field calcu- 
lation using tight binding model, which also takes into 
account the suppression on the d-density-wave order pa- 
rameter at the odd parity surface. 



V 

Wi, i+ i = - 



( C L C i+lj + l/2) - (4,j C *+l,3-l/2) (H) 



into the interaction term (W. The minus-sign between 
the two expectation values takes care of the sign change 
of the d £C 2„ l/ 2-symmetric order parameter under n/2- 
rotation of the lattice. Furthermore for the d x 2_ y 2- 
density-wave state (|10|) we expect W-^+i to alternate 
oc ( — l) 1 , but we do not enforce this behavior in the nu- 
merical mean-field treatment. Using this notation, the 
maximal gap magnitude in the homogeneous bulk state 
is2|W M+ i|. 

Next we exploit the translational invariance along y- 
direction by introducing the partial Fourier transform 
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III. MODEL AND SELF-CONSISTENT 
MEAN-FIELD SCHEME 

As a model for the [110] surface, we consider a strip 
of a rotated square lattice with infinite extension along 
the surface (see Fig. ||). For simplicity we neglect the 
effects of the d-density-wave order on the spin degrees 
of freedom. Then for each spin orientation the kinetic 
energy through hopping between neighboring sites gives 
rise to a term 



Hi 



-t 



(.id) 



K 



c i+l.j + l/2 



J 
-i,3- 



+ c , i c i+l,j-l/2 



(8) 



The index i labels the lattice sites along a zig-zag chain 
in x-direction, while the indices i + l,j ± 1/2 label the 
right upper or lower neighbors of the site i,j. Next we 
introduce a repulsive nearest neighbor interaction 



H V = VJ2 [ c \,^A+x,i+i/2 c 



1,3 



'i+l.j + 1/2 



+ C i ,j Ci ^ C i+l j-l/2 C i+l,J-l/2 



(9) 



We have already left out the interaction between densi- 
ties of spin-up and spin-down electrons on neighboring 
sites, as this does not contribute to the d-density-wave 
decoupling. In real space, the static d-density-wave am- 
plitude between nearest neighbors x and y has the fol- 
lowing form0: 



W(x,y) 



= {4, s ^') = 2 Wo6s ' 8 ' eiQ£ 



I Vy,x — a T" dy,x-\-a 



y,x — b 



W) ( 10 ) 



where d and b denote vectors along the crystal axes with 
length of the lattice constant a. In our rotated system 
(see Fig. |l|), the mean-field decoupling is performed by 
introducing the j- or y-independent imaginary mean-field 



N y denotes the number of lattice sites in y-direction, 
which we choose large enough for the numerical calcu- 
lation. Since the distance in the y-direction between the 
zig-zag chains along the ^-direction is a* = v2a, k = k' 
has to range from —ir/a* to n/a*. In the following we 
set a* = 1. With this the Hamiltonian becomes 



H = - 



>£& 



' 2t > , [ c ik c i+hk + c- + i, fc ca 

i,k 



k 

cos — 
2 



■i .k 
i , k 



C i,k C i+l,k + c i+i : fcCi,fe 



sin ■ 



(12) 



This Hamiltonian can be easily diagonalized for each k 
on the finite strip. The self-consistent scheme is comple- 
mented with the gap equation 



Wi, i+ i = — J2(c\ k c i+ ^ k ) sin - . 



(13) 



From this equation and the eigenfunctions of (H3) we can 
calculate the local d-density-wave amplitude Wa+i- In- 
serting this back into (|12|) we can repeat the process until 
self-consistency is reached. 



IV. NUMERICAL RESULTS 

A. Results for perfect nesting: zero energy states 

Let us first discuss the results at half-filling. We choose 
H = and V = 2t. In the left plot Fig. g we show the spa- 
tial variation of the d-density amplitude at the surface. 
One clearly observes the staggered real space structure 
of the mean-field Wi,t+i, in agreement with dlQ). At 
the surface the amplitude is slightly suppressed on the 
length scale of a few lattice constants, similar to the sup- 
pression of the d-wave superconducting gap function in 



analogous superconducting casesliirc^. In the right plot 
we show the local density of states at the surface and in 
the bulk of the system. The bulk density of states ex- 
hibits a well formed gap edge at 2 1 W^S | with some low 
energy states due to the nodes of the d-density-wave gap. 
The surface density of states shows a pronounced zero en- 
ergy peak. Studying the wave functions corresponding to 
these low energy eigenvalues, we find that the peak in- 
deed arises due to bound states (see lower plot of Fig. ||), 
confirming the expectations from the analysis of Sec. II. 
The bound state wave functions take nonzero values only 
on one sublattice and their decay length depends on the 
size of the d-density-wave bulk gap: for small k' y close 
to the nodes of the bulk d-density-wave gap the bound 
state wave function extends far into the system, while for 
larger values of k' corresponding to a larger bulk gap, 
the decay length is shorter. Furthermore we find surface 
bound states either for positive or negative k' , depending 
on the sign of the d-density-wave amplitude at the sur- 
face. A detailed analysis shows that a sizable part of the 
weight in the zero energy peak is taken from relatively 
high energies « 3£-4£. This appears to be comprehensible 
in view of the small extension of the bound state. 




FIG. 2. Upper left plot: Spatial variation of the (purely 
imaginary) d-density-wave amplitude Wi,i+i. Upper right 
plot: Surface (on the first lattice site, solid line) and bulk 
(dashed line) density of states versus energy. Lower left plot: 
Wave functions of two zero energy bound states correspond- 
ing to transverse wave vectors k' y ~ — 37r/4 (solid line) and 
— 7r/4 (dashed line). The data shown are calculated for a 
system with 300 sites in a:-direction, V = It, /i — and 
T — O.Oli. Lower right plot: Spectrum as function of the 
transverse wavevector k' y . The zero energy states for positive 
and negative k' y are localized on opposite surfaces of the strip. 
All states at nonzero energy extend over the whole system. 



B. Surface bound states, domain walls and domain 
wall bound states away from half filling 

Next we change the band filling away from half-filling 
by varying the chemical potential \x. Then the first strik- 
ing observation is that beyond a certain critical \i away 
from fj, — the self-consistent scheme converges to non- 
homogeneous solutions where the d-density-wave pattern 
exhibits domain walls with phase jumps of 7r, as shown 
in Fig. |3| for /i = — O.li corresponding to an average oc- 
cupation number of ~ 0.98 particles per site. Analyzing 
the charge distribution (rii) (the dashed lines in Fig. || 
show the local deviation from half-filling) , we find that 
the system chooses to form locally half-filled regions with 
intact d-density-wave order. The holes are concentrated 
in the domain walls. With increasing temperature, the 
correlation length of the d-density-wave order gets longer 
and the domain walls evaporate such that the bulk state 
becomes homogeneous again. If we move further away 
from half-filling, the domain wall density increases and 
finally, at another critical chemical potential (determined 
by the strength of the nearest neighbor interaction V, for 
V = 24 at fi « — 0.12i), the local d-density-wave order 
disappears completely. 

No attempt was undertaken to investigate the role of 
the boundary conditions for the domain wall formation or 
to show that these inhomogeneous solutions indeed corre- 
spond to global minima of the total free energy. Further, 
if this was the case, the model should be extended to in- 
clude possible inhomogeneities in the direction along the 
surface as well. Here we do not want to take this difficult 
route, or explore the mean-field phase diagram of this 
model in detail. Instead we will now assume that the do- 
main walls along the [110] direction are stable or can be 
stabilized by some additional interactions, and continue 
to analyze the local excitation spectrum at the surface 
and in the domain wall. 
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FIG. 3. Left plot: Spatial variation of the modulus of 
the d-density-wave amplitude Wi^+i (solid line) and the 
local charge {rii, s ) (dashed line) per spin with respect to 
half-filling.) Right plot: detailed structure of the domain wall 
showing the phase shift by n. The data are calculated on a 
system with 200 sites in a;-direction, V — 2£, /i = — O.li and 
T = 0.01£. 



According to the analytic calculation in Sec. II we 
expect the zero energy bound state to move to finite en- 



ergies. This is also the result of the self-consistent cal- 
culation, as shown in Fig. |j for /i = — 0.08t. For this 
chemical potential the bulk state is still homogeneous. 
Apart from the bound state energy, the whole spectrum 
has moved and is now centered axpund a finita energy 
~ — fj,, as also found by Zhu et al.Ea and WangE3 for the 
impurity case. For the bulk system this shift is equiva- 
lent to the assumption of Chakravarty et al.EJ that upon 
hole-doping away from half-filling the d-density-wave dis- 
persion stays rigid like in a semiconductor and that the 
doped holes form Fermi surface pockets around the nodal 
points of the half-filled state. 

We now include a finite next nearest neighbor hopping 
t'. This has the effect of destroying the (tt,tt) nesting 
in the homogeneous state even at half filling. Neverthe- 
less, we find stable d-density-wave solutions, in particular 
near the so-called van Hove filling for /i = it' where the 
bulk Fermi surface touches the van Hove singularities in 
the density of states at (it, 0) and (0, if). For V = 2t and 
t' = —0.2t, stable solutions exist between fi « — 0.61i 
((n) w 0.92 per site) and fi « -0.85i ((n) « 0.82 per 
site). For band fillings smaller or larger than these criti- 
cal values the mean-field solutions for the d-density-wave 
state collapse rapidly with very narrow doping regions 
close to these critical dopings where inhomogeneous solu- 
tions occur. This suggests that the nonzero next nearest 
neighbor hopping t' inhibits the formation of well shaped 
domain walls with regions of intact d-density wave order 
and charge accumulations in between. In the right plot 
of Fig. [| we plot the surface and bulk density of states 
for the homogeneous case t' = —0.2t and fj, = —Q.7t. 
Both curves have a large asymmetry with respect to zero 
energy and the peak in the surface density of states is 
centered at comparably high excitation energies. 






Energy / 1 



-0.5 0.5 

Energy / 1 



FIG. 4. Left plot: Local density of states at the surface 
(solid lines) and in the bulk density of states (dashed lines) 
versus energy, for V = It, \x = -0.08* and T = O.Oli 
(thin lines) and T = 0.08i (thick lines), normalized with re- 
spect to the values above the d-density wave transition at 
T c w m 0.125t. The local density of states was smeared out 
by Lorentzians with width 2T. Right plot, upper curves: sur- 
face (solid line) and bulk (dashed line) density of states in 
the d-density-wave state for t' = —0.2t and fi = — 0.7t at 
T = 0.02i. Right plot, lower curves: surface (solid line) and 
bulk (dashed line) density of states for homogeneous d-wave 
superconducting order coexisting with d-density-wave order 
(non-self-consistent) of comparable amplitude. 



Now let us turn back to the domain walls for t' = 
at sufficient doping away from half filling. As shown in 
Fig. ||, the local density of states at the domain wall 
exhibits a low energy peak similar to the one at the sur- 
face. For a given transverse wavevector we can analyze 
the wave functions contributing to this peak. Two of 
these wave functions for different transverse wavevectors 
k' are shown in the right plot of Fig. |5|: the lower curve 
belongs to a wave function with k' a* = —3, while the 
upper curve is for k' a* = —2.5. The latter wave function 
has visible amplitude in two domain walls. The difference 
between the two wave functions can be easily understood: 
the wave function with larger transverse momentum be- 
longs to a larger bulk gap (in our coordinate system, the 
bulk gap is maximal for k' a* = n). Therefore the wave 
function decays rapidly away from the domain wall. The 
wave function for the smaller k' has a larger decay length 
and thus can effectively hybridize with solutions localized 
at the neighboring domain wall. For even smaller values 
of the transverse wavevector, the wave functions corre- 
sponding to the sub-gap energy states distribute their 
weight over several domain walls. 
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FIG. 5. Left plot: Local density of states at a domain wall 
(solid line) and bulk density of states (dashed line) versus 
energy The data are calculated on a system with 300 sites in 
^-direction, V - 2t, H = -0.lt and T = O.OOli. The local 
density of states was smeared out through Lorentzians with 
width O.Olt. Right plot: The upper (k' y a* = —2.5) and lower 
(k' y a* = —3) curve (shifted by ±0.5) show wave functions of 
states with different transverse wavevector k contributing to 
the low energy peak at a domain wall. The thin line in the 
middle shows the absolute magnitude of the d-density-wave 
amplitude Wi,i+i. 



V. CONCLUSIONS 

We have calculated the local density of states in the 
presence of d-density wave order on the two-dimensional 
square lattice. At [110] surfaces we find zero energy 
(for perfect (it, n) nesting of the Fermi surface at half- 
filling) or low energy (away from perfect nesting) sur- 
face states, which give rise to peaks in the calculated 
surface density of states, in good qualitative agreement 
with recent.res.ults for the local density of states around 
impuritiesEflEl The surface density of states is mea- 
surable in in-plane tunneling experiments, as demon- 
strated for the high-T c cuprates in the superconduct- 



ing stat< 



Similar to the proposals by Zhu et al 



, , , , -.fi 

WangE3 and MorrEj, this provides another experimen- 
tal test for possible d-density-wave long range order in 
the pseudogap phase of the underdoped cuprates, as pro- 
posed by Chakravarty et alJj. Since surface effects cru- 
cially depend on the quality of the surface, it could still 
be that the bound state features are wiped out by surface 
imperfections, even if d-density-wave order is present. 
Moreover the temperature smearing reduces the low bias 
peaks considerably. Therefore we propose to perform 
experiments on strongly underdoped samples with large 
pseudogap and low superconducting transition tempera- 
tures which have [110] surfaces that exhibit a zero-bias 
peak in the superconducting state. The effect of surface 
imperfections on the latter feature should be comparable 
to that on the d-density bound states. Thus, if the su- 
perconducting state does exhibit an Andreev bound state 
one may expect a d-density surface state in the pseudo- 
gap phase as well. We note that so far the experimental 
results do not point towards d-density wave order: pub- 
lished conductance measurements on underdoped YBCO 
films and LSCO single crystalaU and tunneling .data on 
001 tilt grain boundaries in underdoped YBC013 show 
zero energy peaks below the superconducting T c and no 
feature above T c . 

For finite next-nearest neighbor hopping, the peak in 
the calculated local density of states due to the surface 
states moves to rather large excitation energies, such that 
one might wonder if this peak could still be distinguished 
from a gap edge similar to the one in the bulk density 
of states. We think that the distinction is still possible, 
because the peak due to the surface states has to move 
to zero energy when the system enters the <i-wave su- 
perconducting state. In contrast with that the gap edge 
feature due to the extended bulk states should move to 
higher excitation energies in the superconducting state 
(see lower curves in the right plot of Fig. |j), if the latter 
is indeed formed on top of the d-density-wave order, as 
suggested by Chakravarty et al.lil. One may even try to 
observe a shift of the bound state feature - if existent 
above the superconducting T c - from nonzero energies 
due to the imperfect nesting to zero energy when one 
goes from the d-density-wave phase into the d-wave su- 
perconducting state. Another interesting experimental 
question might be the search for domain walls of the d- 
density-wave order parameter along the [110] direction 
which again should give rise to low energy peaks in the 
local c-axis tunneling density of states. We note however 
that our calculation is not a stringent proof of the exis- 
tence of domain walls for the reasons explained above. 

It is a pleasure to thank Chetan Nayak, Ying Liu, Masat- 
sura Igami, Manfred Salmhofer and Katsunori Wak- 
abayashi for helpful discussions. This work was finan- 
cially supported by the Swiss National Science Founda- 
tion. 



S. Chakravarty, R.B. Laughlin, D.K. Morr, and C. Nayak, 
Phys. Rev. B 6 3, 094503; S. Chak ravarty, Hae- Young Kee, 



2 

3 C 

4 



and C. Nayak 
H.J. Schulz, Phys 
Nayak, Phys 



ond-mat/0101204 



Rev. B 39, 2940 (1989). 
Rev. B 62, 4880 (2000). 
I. Affleck and J. B. Marston, Phys. Rev. B 37, 3774 (1988); 
I. Affleck, Z. Zou, T. Hsu, and P.W. Anderson, Phys. Rev. 
B 38, 745 (1988). 

5 X.G. Wen and P.A. Lee, Phys. Rev. Lett. 76, 503 (1996); 
J. Kishine, P.A. Lee, and X.G. Wen, Phys. Rev. Lett. 86, 
5365 (2001). 

6 J. O. Fjaerestad, J. B. Marston, |cond-mat/0107094 

7 Qiang-Hua Wang. Jung Hoon Han, and Dung-Hai Lee, 
cond-mat/0011398 . 

8 Wonkee Kim, Jian- Xin Zhu, J. P. Carbotte, and C. S. Ting, 



cond-mat/0106594 



u D.J. Van Harlingen, Rev. Mod. Phys. 67, 515 (1995). 

10 M. Sigrist, Prog. Theo. Phys. 99, 899 (1998). 

11 For a review, see S. Kashiwaya and Y. Tanaka, Rep. Prog. 
Phys. 63, 1641 (2000). 

12 M.Covington, M.Aprili, E.Paraonu, L.H.Greene, F.Xu, 
J. Zhu, and C.A. Mirkin, Phys.Rev.Lett. 79, 277 (1997). 

13 L. Alff, S. Kleefisch, U. Schoop, M. Zittartz, T. Kemen, 
T. Bauch, A. Marx, and R. Gross, Eur. Phys. J. B 5, 423 
(1998). 

14 Y. Dagan, A. Kohen, A. Revcolevschi, and G. Deutscher, 
Phys. Rev. B 61, 7012 (2000). 

15 N.C. Yeh, J.Y.T. Wei, C.T. Chen, W.D. Si, and X.X. Xi, 
Physica C 341, 1639 (2000). 

16 C. Walti, H.R. Ott, Z. Fisk, and J.L. Smith, Phys. Rev. 
Lett 85, 5258 (2000). 

17 F. Laube, G. Goll, H. v.Lohneysen, M. Fogelstrom, and F. 
Lichtenberg, Phys. Rev. Lett. 84, 1595 (2000). 

18 Jian-Xin Zhu, Wonkee Kim, C.S. Ting, and J. P. Carbotte, 



cond-mat/010558 



Qiang-Hua Wang 
D.K 



Morr, 



I 



oiKl-ma.t/0106327| . 



ond-mat/0106285 



Adv. 



See for example M.S. Dresselhaus and G. Dresselhaus, 
Phys. 30, 139 (1981). 

A similar approach to graphite edge states has been dis- 
cussed by M. Igami, unpublished. 

M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusakabe, 
J. Phys. Soc. Jpn. 65, 1920 (1996). 

C. Honerkamp, K. Wakabayashi, and M. Sigrist, Euro. 
Phys. Lett. 50, 368 (2000). 



